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Abstract 
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I. INTRODUCTION 



When a particle is exposed to a strong potential field, the relativistic effect must be 
considered which gives the correction for nonrelativistic quantum mechanics. Taking the 
relativistic effects into account, a spinless particle in a potential field should be described 
with the Klein-Gordon (KG) equation. The solution of the Dirac equation is considered 
important in different fields of physics like nuclear and molecular physics [1,2]. Within 
the framework of the Dirac equation the spin symmetry arises if the magnitude of the 
spherical attractive scalar potential S and repulsive vector V potential are nearly equal such 
that S ~ V in the nuclei (i.e., when the difference potential A = V — S = C- = Ca, 
with Ca is an arbitrary constant). However, the pseudospin (pspin) symmetry occurs if 
S ~ —V are nearly equal (i.e., when the sum potential E = V + S = C + = Cs, with 
Cs is an arbitrary constant) [3]. The spin symmetry is relevant for mesons [4]. The pspin 
symmetry concept has been applied to many systems in nuclear physics and related areas 
[3-7] and used to explain features of deformed nuclei [8], the super-deformation [9] and 
to establish an effective nuclear shell-model scheme [5,6,10]. Recently, the spin and pspin 
symmetries have been widely applied on several physical potentials by many authors [11-26]. 
For example, the Dirac equation has been solved for the deformed generalized Poschl- Teller 
(PT) potential [27], modified PT potential [28,29], Manning-Rosen (MR) potential [30], well 
potential [31], modified Rosen- Morse (RM) potential [32] and class of potentials including 
harmonic oscillator, Morse, Hulthen, Scarf Eckart, MR, Trigonometric RM potentials and 
others [33] in the framework of the approximation to the spin-orbit centrifugal term using 
the proper quantization rule, algebraic methods, Ladder operators and su(2) algebra. 

The exact solutions of the Dirac equation for the exponential-type potentials are possible 
only for the s-wave (k = ±1 case) when the spin-orbit coupling term will get suppressed 
[34]. However, for /-states an approximation scheme has to be used to deal with the spin- 
orbit centrifugal k(k+ l)/r 2 (pseudo-centrifugal, k(k — l)/r 2 ) term. In this direction, many 
works have been done to solve the Dirac equation with large number of potentials to obtain 
the energy equation and the two-component spinor wave functions [35-42]. It has been 
concluded that the values of energy spectra may not depend on the spinor structure of the 
particle [43], i.e., whether one has a spin-1/2 or a spin-0 particle. Also, a spin-1/2 or a spin-0 
particle with the same mass and subject to the same scalar S and vector V potentials of 
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equal magnitude, i.e., S = ±V (A = E = or C± = 0), will have the same energy spectrum 
(isospectrality) , including both bound and scattering states [43]. It has been shown that for 
massless particles (or ultrarelativistic particles) the spin- and pspin spectra of Dirac particles 
are the same for the harmonic oscillator potentials [44]. 

Recently, we obtained the spin symmetric and pspin bound state solutions of the Dirac 
equation with the standard RM well potential model [20,45]: 

V(r) = —V\ sec h 2 ar + V2 tanh or, (1) 

where the coupling constants V\ and V2 denote the depth of the potential and a is the range 
of the potential that has an inverse of length dimension. We use the computer software 
MATLAB and plot the potential (1) for three different set of parameters V\ and Vi- It is 
plotted in Fig. 1. 

The aim of the present paper is to extend the s-wave solutions by solving the Dirac equa- 
tion with some physical potentials given in Ref. [34] in the framework of the Nikiforov-Uvarov 
(NU) method [46-50] by taking an approximation to deal with the centrifugal (pseudo- 
centrifugal) potential term [20,51]. The approximation scheme used to deal with the spin- 
orbit centrifugal barrier k(k + l)/r 2 holds for values of spin-orbit coupling quantum number 
k that are not large and vibrations of the small amplitude [51]. In the presence of spin sym- 
metry S ~ V and pspin symmetry S ~ —V, we calculate bound state energy eigenvalues 
and their corresponding upper and lower spinor wave functions. We also show that the spin 
and pspin symmetry Dirac solutions when A = C_ and E = C + can be reduced to the exact 
spin symmetry and pspin symmetry limitation A = and E = 0, respectively. These are 
found to be identical to the KG solution for the V = ±S cases. Furthermore, the bound 
state solutions of the Schrodinger equation are also obtained from the nonrelativistic limit 
of the Dirac equation if an appropriate mapping of parameters is used. 

The paper is organized as follows: Section 2 is mainly devoted to the basic spin and 
pspin Dirac equation. In sect. 3, the approximate analytical bound state solutions of the 
(3 + l)-dimensional Dirac equation with the reflectionless-type and the RM potentials are 
obtained in the presence of the spin and pspin limits using a parametric generalization of the 
NU method. In sect. 4, special cases like the s-wave k — ±1 (1 = 1 = 0) and nonrelativistic 
limit are studied. Section 5 gives the relevant conclusion. 
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II. BASIC SPIN AND PSPIN DIRAC EQUATIONS 



The Dirac equation for fermionic massive spin-1/2 particles subject to vector and scalar 
potentials is [1] 

[ca ■ p + f3 (Mc 2 + S(r)) + V(r) - E] ^ nK (v) = 0, ^(r) = ^(r, 9, 0), (2) 

where E is the binding relativistic energy of the system, M is the mass of a particle, p = 
—iftV is the momentum operator, and a and (5 are 4x4 Dirac matrices [3,12,20]. The spinor 
wave functions take the form 

. M l( Fn.{r)Y> m {9A) \ _ 
Vw(r) = - , 3 

r \iG nK {r)Y l jm {9A) ) 

where F nK (r) and G nK (r) are the radial wave functions of the upper- and lower-spinor compo- 
nents, respectively, and Y- m (9, <fi) and Yj m (9, <f>) are the spherical harmonic functions coupled 
to the total angular momentum j and it's projection m on the z axis. 

In the presence of spin symmetry ( i.e., A = C_ = Ca), one obtains a second-order 
differential equation for the upper-spinor component [12,20,52,53]: 

KJr) ~ + Al + B s ^ F nK (r) = 0, (4) 

where 

Al = ^ [MV - El - (Mc 2 - E nK ) C_] ,B s = ^- 2 (Mc 2 + E nK - C_) , (5) 

and k(k + 1) = 1(1 + 1), k = I for k < and k — — (I + 1) for k > 0. Further, the 
lower-spinor component can be obtained as 

g -^= M* + L-cM + ^) F "« {r) ' (6) 

where E nK ^ —Mc 2 when C_ = Ca = (exact spin symmetric case). It means that only 
positive energy spectrum is permitted. 

Overmore, in the presence of pspin symmetry ( i.e., S = C + = Cs), one obtains a 
second-order differential equation for the lower-spinor component, 

G'L(r) - + A 2 - B ps a] G nK (r) = 0, (7) 



4 



where 



A 2 



1 



[M 2 c 4 - E 2 nK + (Mc 2 + E nR ) C + ] , B p 



1 



h 2 c 2 L- -nn ■ V — ■ , ^2 

The upper-spinor component F nK (r) can be obtained by means of 

1 



(Mc 2 - £ nK + C+) . (8) 



(9) 



Mc 2 - £ 

where i? nre 7^ Mc 2 when C + = Cs = (exact pspin symmetric case). It means that only 
negative energy spectrum is allowed for this case. From the above equations, the energy 
eigenvalues depend on the quantum numbers n and k, and also the pseudo-orbital angular 
quantum number I according to k(k — 1) = 1(1 + 1), which implies that j — I ± 1/2 are 
degenerate for / 7^ 0. The quantum condition for bound states demands the finiteness of the 
solution at infinity and at the origin points. 

It is known that Eqs. (4) and (7) can be solved exactly only for the case of k — — 1 
(I = 0) and k — 1 (I — 0), respectively when the spin-orbit coupling centrifugal and pseudo- 
centrifugal terms will get suppressed. In the case of nonzero I or I values, we can use the 
approximation scheme to deal with the spin-orbit centrifugal (pseudo-centrifugal) term when 
k is not large and when vibrations of the small amplitude near the minimum point r = r e 
[17,51] 



D + D 1 - eXPi 7 2 ^\ +D 2 



— exp(— 2ar) \ 2 
1 + exp(— 2ar) J 



1 + exp(— 2ar) 

where A is the parameter of coefficients (i = 1,2,3) given by 

A. = l-( : T \ , x " (3 + 2«r e ) 



D x = -2 (exp(2ar e ) + 1) 



D 2 = (exp(2ar e ) + 1) 



2ar e J \ 1 + exp(— 2ar e 
1 +exp(-2ar ey , 

(3 + 2ar e ) 



2ar e 
1 + exp(— 2ar e 



1 + exp(— 2ar e ) 
2ar P 



2ar P 



3 + 2ar e - 



Aar e 
1 + exp(— 2ar e ) 



(10) 

(11a) 
(lib) 
(11c) 



and higher order terms are neglected. 



A. Spin symmetric solution 



We take the sum potential in Eq. (4) in the form of standard RM well potential (1), i.e., 



E = V(r) = -4Vl 



exp(— 2ar) 
(1 + exp(-2ar)f 



+ V 2 



1 — exp(— 2ar)\ 
1 + exp(— 2ar) ) 



(12) 
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Upon introducing the new variable z(r) = exp(— 2ar) and substituting the above sum 
potential into Eq. (4) which then can be cast into the form 

Z{l + Z) z z (1 + z) 

where at boundaries we require that F nK (0) = F nK (—l) = and the parameters (i — 0, 1, 2) 
take the forms: 



a =' ^-^D + B s V 2 + A 2 >0, 



e 



4a 2 V r 



e 



1 /k(k + 1) 



a 2 = ^ ^^7^ (A) - + A) + ^ - ^ J • (14) 
We begin the application of the NU method [46-50] by comparing Eq. (13) with the hyper- 



geometric differential equation 

€(r) + ^<(r) + ^Mr) = 0, (15) 

where 

*Pn(r) = (f)(r)y n (r), (16) 

to identify the parameters, 

t(z) — 1 + z, <j(z) — z (1 + z) , 5(2) = — a 2 z 2 + a\z — a , (17) 
and further calculate the function n(z) as 



lr ... „ . /I r . M2 



= 2 [ CT '( r ) - r ( r )l ± y 4 K( r ) - r ( r )l - a ( r ) + ka ( r ) 



= Z - ± i J[l + 4 (a 2 + k)] z 2 + 4 (fc - ai) z + 4a 2 ,. (18) 

Now we also seek for a physical value of /c that makes the discriminant of the expression 
under square root in Eq. (18) to be zero, that is 

k — di + 2oq ± 2a q, 



+ + (19) 
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Upon the substitution of the value of k into Eq. (18), we obtain the following convenient 
solutions: 

(20) 



tt(z) = -oq - - (2a + q-l)z, 



and 



k — a± + 2oq + a q. (21) 

With regard to Eqs. (17) and (20), we can calculate the function t(z) = t(z) + 2tt(z), taking 
into consideration the bound state condition which has to be established when t'(z) < 0, as 



t(z) = 1 - 2a - (2a + q-2)z, r'(z) = - (2a + q - 2) < 0, 



(22) 



where prime denotes the derivative with respect to z. According to the method [46-50], in or- 
der to find the energy equation from which one calculates the energy eigenvalues, we need to 
find the values of the parameters: A = k+n'(s) and A — X n — — nr'(s)- \n (n — 1) cr"(s), n = 
0, 1, 2, • • • , as 

1 / 1 \ 

(23) 



A = ^ + ai + 2a 2 Q - a + ( a - - J <j. 



and 



X n = -n 2 -n + n (2a + q) , n = 0, 1, 2, • ■ • , (24) 

respectively. Using the relation A = A„ and the definitions of variables in Eqs. (14) and 
(19), we obtain the transcendental energy equation of relativistic spin- 1/2 particles in the 
presence of vector and scalar potential, 

A K^ 4 - EL - (Mc* - E nK ) C_] = - ^t^ " B S V 2 



h 2 c 



+ a z 



1 q K{K+l){D 1 -D 2 )/rl + 2B s V 2 
U+ 2 2 + 4a=»(n + I-§) 



-I 2 



(25) 



Furthermore, in the exact spin symmetric case (i.e., V = S, A = 0, C_ — > 0), we obtain the 
arbirary /-wave energy equation in the KG theory with equally mixed RM-type potentials 
(in units of h — c — 1), 



M2 _ ^ = JV+gDo 



(E nl + M) V 2 



+ a z 



n + 5 + 



1(1 + 1) (A - D 2 ) jrj + 2 (E nl + M) V 2 
Aa 2 (n + 5) 



(26) 
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with 



err* cr 



(27) 



where the quantum number n — 0, 1, 2, • • ■ , and the orbital quantum number / = 0, 1, 2, • • • . 
Actually, the above expression resembles Eq. (13) reported in [54] when I = (s-wave case). 
Now, we are going to find the corresponding wave functions for the present potential model. 
Firstly, we calculate the weight function defined as 



and the first part of the wave function in Eq. (16) as 

(j>{z) = exp ^ ^d^j = z ~ a ° (! + z Y l[1 ~ q) ■ 



(28) 



(29) 

Hence, the second part of the wave function in relation (16) can be obtained by means of 
the so called Rodrigues representation 

yniz) = K(r)p(r)] = Knz2ao (1 + z)q ^ [zn ~ 2ao (1 + z)n ~ 9] 

„ P i-**,-*)(l + 2z), ze [0,1], (30) 

where the Jacobi polynomials Pi fJ, ' u \x) are defined for Re(v) > —1 and 3?(/i) > —1 for the 
argument x G [—1, +1] and K n is the normalization constant. By using F nK (z) = (f)(z)y n (z), 
in this way we may write the upper-spinor wave function in the fashion 

F nK (r) = K nK (exp(-2ar)r ao (1 + exp^ar))^ 1 ^ pH** -«)(1 + 2exp(-2ar)) 

= N nK (exp(-2ar))~ a ° (1 + exp(-2ar))5 (1_9) 2 F X (-n, n + 1 - 2a - q; -2a + 1; exp(-2ar)) , 

(31) 

where a > 0, q > — 1. The calculated normalization constants K nK for the upper-spinor 
component are 

-1/2 

, (32) 



N = 



T(-q + 2)r(-2a + 1) ^ (-l) m (n - 2a + 1 - q) m T (n + m) 



2«r(n) 



m=0 



ml (m — 2a )!T (m — 2a — q + 2) 



with 



fun = 3 F 2 (-2a + m, -n, n + 1 - 2a - q; m - 2a - q + 2; 1 - 2a ; 1) . 



(33) 
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In addition, the corresponding lower component G nK (r) can be obtained as follows 



(exp(-2ar))~ a ° (1 + exp^ar))^ 1 "') 



a (1 — q) exp(— 2ar) k 

2aa ° — rr> 1 — o — ^Y~~ + ~ 

(1 + exp(— 2ar)) r 



+c r , 



(Mc 2 + E nK - C_) 
x 2 Fi {—n, n — a — q + 1; — 2a + 1; exp(— 2ar)) 
2an [n-2a -q + 1] (exp(-2ar))" ao+1 (1 + exp(-2ar))^ (1 ~ g) 



(2a + 1) (Mc 2 + E nK - C_) 

x 2 Fi (-n + 1; n - 2a - g + 2; -2a + 2; exp(-2ar)) , a > 0, (34) 

where E nK ^ —Mc 2 for exact spin symmetry. Here, note that the hypergeometric series 
2 Fi (—n, n — 2a — q + 1; — 2a + 1; exp(— 2ar)) terminates for n = and thus converge for 
all values of real parameters q > and a > 0. 



B. Pspin symmetric solution 

In the same way as before, this time taking the difference potential in Eq. (7) as 

A = V(r) = -4K ex P(~ 2 ^) 2 + y 2 f l-«P(-2«rA (35) 
V; (l + exp(-2w)) 2 V 1 + exp(-2ar) J V ; 

and in terms of new variable z(r) = exp(— 2ar), leads us to obtain a Schrodinger-like equa- 
tion for the lower-spinor component G nK (r), 

Z{L + Z) z l (1 + z) 

where the parameters bj (j = 0, 1, 2) are defined by 



bo = ^J^^D + B ps V 2 + Al s >0, 



e 



h = -L (A - 2 D ) + 43^ - 2^. 

& 2 = ^ (A> - A + I> 2 ) + < - B ps V 2 ^j . (37) 

To avoid repetition in the solution of Eq. (36), a first inspection for the relationship between 
the present set of parameters (&o,&i,&2) and the previous set (00,01,02) tells us that the 
negative energy solution for pseudospin symmetry such that £ = C + = Cs can be obtained 
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directly from those of the positive energy solution above for spin symmetry by performing 
the changes [12,20]: 

F nK (r) o G nK (r), V(r) -+ -V(r) (or V 1 -> -V u V 2 -> -V 2 ), E nK -> -E nK and C_ -C 4 

(38) 

Considering the previous results in Eq. (25) and applying the above transformations, we 
finally arrive at the pspin symmetric energy equation 



[M 2 c 4 - El K + (Mc 2 + E nK ) C + ] = 



H 2 c 2 k(k-1)D 



+ (Mc 2 - E nK + C+) V 2 



h 2 c 2 a 2 



+ 



where 



2n + 1 - p + 



h 2 c 2 K (k - 1) (A - £ 2 ) /r e 2 - 2 (Mc 2 - £ nK + C+) V 2 
h 2 c 2 a 2 a 2 (2n + 1 — p) 



P 



i+ n{K-Y)D 2 _W } Bp 1 



a 2 r 2 



a" 



, (39) 



(40) 



Again, the radial lower-spinor wave function in Eq. (31) becomes 

G nK (r) = d nK (exp(-2ar)r b0 (1 + exp(-2ar))* (1 - p) P^ 2b ^{l + 2 exp(-2ar)). 

= d nK (exp(-2ar))~ 6 ° (1 + exp(-2ar))^ (1 " p) 2 F X (-n, n + -2b -p+1; -2b + 1; exp(-2ar 

(41) 

which satisfies the restriction condition for the bound states, i.e., p > and b > and the 
normalization constants is 



d n 
with 



T(-p + 2)r(-2fe + 1) ^ (-l) m (n + 1 - 26p - p)) m I> + m) 



-1/2 



, (42) 



2ctT(n) ^ m!(m-26 )!r(m-26 -p + 2) 

3™ = 3 F 2 (-26o + m, -n, n + 1 - 26 - p); m - 2b - p + 2; 1 - 26 ; 1) . (43) 



III. APPLICATIONS TO SOME PHYSICAL POTENTIAL MODELS 

We adopt the following two physical potential cases that belong to the general potential 
model been introduced in Eq. (1). 
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A. The reflect ionless-type potential 



The reflectionless-type potential is the special case of the symmetrical double-well poten- 
tial. offered by Biiyiikkihg et al [55] to describe the vibration of polyatomic molecules. This 
can be achieved when the coefficient of tanh ar becomes zero. So, it takes the form 

V(r) = -a 2 sec h 2 ar, a 2 = A(A + l)/2, A = 1, 2, 3, • • • . (44) 

The potential is plotted in Fig. 2 for three different values A = 1,2 and 3. It follows that 
the energy equation in Eq. (25) becomes 



2 



M>-El K -CAM-E m ) = ->±^ + ^ 



2n + 1 - g + 



k (k + 1) {D 1 - D 2 ) /r 



a 2 (2n + 1 - g ) 



5) 



and the upper-spinor wave functions from Eq. (31) turns to be 

F nK (r) = N nK (exp(-2ar))- so (1 + exp(-2ar))^ (1 -' ?o) pM^) (1 + 2 exp(-2ar)), (46) 
where 



1 / «(« + !) „ /, , k(k + 1)D 2 , 4a*B a 

2^\M^ D0 + As> °> % = \ l+ «V? + — 



e 



It is worth noting that the results given above in Eq. (45) and (46) are identical to those 
ones of Ref. [27] for s-wave case (k = —1). In the presence of the pspin case, the energy 
spectrum becomes 



2 L -2 r, l' \ k(k-1)D , c? [ o _ 1 _ , K(K-l)(D 1 -D 2 )/rl 

( 



M 2 -El + C + (M + E nK ) = - ^ ;^° +^ 



2n + 1 - p + 



a 2 (2n + 1 - p ) 



8) 



r 2 4 

and the lower spinor component of pseudospin symmetric wave function 

G nK (r) = N nK (exp(-2ar))- wo (1 + exp(-2ar))^ (1 ~ po) P^ 2w ^\l + 2 exp(-2ar)). (49) 
where 



1 / «(«-!) _ / «(«-!) A 4a 2 5 ps 



e V e 



4a 2 £? ps /a: 2 < 1 for bound states when k — 1. 

Let us now discuss the non-relativistic limit of the energy eigenvalues and wave functions 
of our solution. If we take C_ = (A = 0) and consider the transformations E nK + M ~ 2/i 

11 



and E nK — M ~ E ni [52,53], we would have the following expression for the energy equation 
(45) and wave functions (46) (in h — c — 1) 

2 



E, 



nl 



1(1 + 1) D a 2 



2/ir 



2fi 



1 1 m (/ + 1) (D l - D 2 ) 

n + o ~ + 777/ — 71 7T 

2 2 4r2« 2 (n + | - ±q ) 



(51) 



and the wave functions: 



Rni(r) = N nl (exp(-2ar))- so (1 + exp(-2ar)p (1 ^ o) ^""'(l + 2 exp(-2ar)), (52) 



with 



2a 



h 2 



a 2 r 2 



(V 



h 2 ' 



(53) 



where E n \ < D is a condition for bound state solutions. 2. 

To conclude, it is necessary to mention that the reflectionless-type potential here reminds 
one of the modified PT potential in the one-dimensional (ID) case [29]. However, for the 
present case it is in the three-dimensional (3D) case. Thus, the original symmetry is broken. 
The energy levels could be obtained readily. 



B. The Rosen-Morse potential 

The standard RM potential was given by Rosen and Morse in Ref. [45] useful to describe 
interatomic interaction of the linear molecules and helpful for discussing polyatomic vibra- 
tional energies. As example of its application to the to the vibrational states of the iVi/3 
molecule. This can be achieved when 



V(r) = —a(a + a) sec h 2 ar + 2b tanh ar, 



(54) 



where a and b are real dimensionless parameters. In Fig. 3, we plot this potential for three 
various sets of parameter values. It follows that from Eqs. (25) and (31), the spin symmetry 
energy spectrum for the RM well is 

M 2 -E 2 nK -C- (M - E nK ) = — — + 1} Do — 2b (M + E nK - C.) 



+ 



cr 



2n + 1 - q 1 + 



k(k + 1) (Pi - D 2 ) /r 2 + 4b (M + E nK - C_) 
a 2 (2n + 1 - qi) 



(55) 
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and the upper spinor component F nK (r) of the wave functions as 

F nK (r) = N nK (exp(-2ar)y Sl (1 + exp(-2ar))^ (1 -' ?l) P^ 2s ^\l + 2 exp(-2ar)), (56) 

respectively, where 

1 /«(« + 1) ~ ~~~~ ~ rt /" k(k + 1)D 2 4a(a + a)B s 

s > = *i J ir LD ° + 2bB ' + A - > °- 91 = V 1 + "W - + -4^- (57) 

Overmore, in the presence of the pspin symmetry, the energy spectrum for the RM well is 



M 2 - £ 2 K + C + (M + £ nK ) 



2&(M-£ nK + C + ) 



+ T 



2n + 1 — pi H 



a 2 (2n+ 1 - pi) 



(58) 



and the lower-spinor wave function is 



G nK (r) = d nK (exp(-2ar))^ 1 (1 + exp(-2ar))* 



i(l-Pi) p(-2tui,-pi) 



(1 + 2exp(-2ar)), (59) 



with 



I k ^_A Dq + 2bB ps + A 2 S > 0, Pl = A / 1 - 



«(k-1)D 2 4a(a + ai)-Bp £ 



a 2 r 2 



cr 



where 4a(a + a)B ps /a 2 < 1 when k = 1. Let us now discuss the non-relativistic limit of 
the energy eigenvalues and wave functions of our solution. If we take C_ = (A = 0) and 
consider the nonrelativistic limits [52,53], we would have the following expression for the 
energy equation (55) and the upper spinor component of the wave functions (56) (in units 

h = c = \) 



2/xr, 



or 
'2/x 



-1 2 



l(l + l)D 2 8fia(a + a) 



1(1 + 1) (D l - D 2 ) jr 2 e + 8 fib 



a 2 r 2 



cr 



a 2 [n+ h ~ k\/l + 



l(l+l)D 2 _|_ 8iia(a+a) 



(61) 



and 



Rni{r) = N nl (exp(-2ar))^ Sl (1 + exp(-2ar)p (1 - ?l) pH 2 * -90(1 + 2 exp(-2ar)), (62) 



respectively, where 



51 = Ya 



1 ^^-D + A/ib - 2/iE nl > 0, qi 



1 + W + W (63) 



a 2 r 2 



cr 



13 



and N n i is the normalization constant. 

To conclude, it is necessary to mention that the RM potential was studied byusing the 
proper quantization rule in Ref. [33]. 



IV. DISCUSSIONS 



We study two special cases of the energy eigenvalues given by Eqs. (25) and (39) for the 
spin and pspin symmetry, respectively. 

(I) The s-wave spin symmetric case (re = — 1, I — 0) (in units h — c — 1). For the 
reflectionless-type potential, we have 



M 2 - E 2 - (M - £ n ,-i)C_ = a 2 



2 2 



I 1 / 4a 2 
" + + — (M + E n -i - C- 



(64) 



and 



|(i-Vi+^(Af+^,-i-C-)) 



F n ,-i(r) = c n ,_! (exp(-2ar))- A - l/2a (1 + exp(-2ar)) : 

(-A^/a-Jl+^-iM+E^-C-)) 

xP n v y (l + 2exp(-2ar)) (65) 

where A 2 _ 1 = M 2 — E 2 _ 1 — (M — E n _i) C_. The s-wave pspin symmetric case (re = 1, / = 0): 



M 2 - < x + C + (M + K,i) = 



1 1 / An 2 



(66) 



where ^ (M — E nA + C + ) < 1 and the lower spinor component of pspin symmetric wave 
function 

G B>1 (r) = d nA (eM-2ar)y A ^ (1 + exp(-2«r))K 1 -^^^) 

(-A 1 /a,-Jl-*g(M-E nl +C+)) 

xP n v y (l + 2exp(-2ar)). (67) 

where A\ = M 2 — E\ x + (M + E n ^) C+. For the RM potential model (spin symmetric case), 
we have 

M 2 - E 2 n ^ — C- (M — K,-i) = -26 (M + P n ,-i - C_) 



+ of 



n 2 



1 b(M + E nK -C. 

n + — + 



2 2 



(»+i-¥)J ' 



(68) 
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and the upper spinor component F nK (r) of the wave functions as 

F n _i(r) = iV n _i (exp(-2ar))~ 7 - 1 / 2 (1 + exp(-2ar))^ (1 " /3 - l) p^-T-i -^-i)(i + 2exp(-2ar)), 



4a(a+«) 
a* 

For the pseudospin case: 



where = Jl + (M + F n ,-i - C-) and 7 -i = v 7 ^ (M + £ n ,_i - C_) + A 2 _ 1 /2a. 



M 2 - P 2 a + C + (M + £ n>1 ) = 2b (M — E nA + C+) 

(69) 



1 6(M-F n;1 + CV ' J 

n H 1 



+ " 2 ' 2 2; ^(n+|-%) 
and the lower-spinor wave functions become 

G n>1 (r) = 4,i (exp(-2ar))~ 7l/2 (1 + exp(-2or))* (1_/Jl) + 2 exp(-2ar)), (70) 

where ft = ^/l - (M - £ n ,i + C+) and 7l = ^2b(M - E n , l + C + ) + A{/a. 

(II) The transformation of the potential (1) into other potential forms. For a potential 
V(x), when one makes the transformations: x — > —x, a — > ia and V2 iV~2 (complex 
parameters), then Eq. (1) transforms into a trigonometric Rosen- Morse-type (tRM) form: 

V(x) = — V\ sec 2 ax + V2 tan ax, a = — , x=[0,a], (71) 

where 5R(Vi) > 0. When x —> — x and % — > — i, if the relation V(— x) = V* exists, the potential 
V(x) is said to be PT-symmetric, where P denotes parity operator (space reflection) and 
T denotes time reversal (see e.g., [56,57] and the references therein). This PT-symmetric 
potential is plotted in Figure 4 for various sets of parameters V\ and V-i- Thus the spin- 
symmetric energy equation (k = —1) can be obtained from Eqs. (19) and (25) as 

M 2 - - (M - £„,_,) C_ = -a 2 \ - ^ (M + £„,_, - C_) 



2a 



where 4Vl (M + P n _x - C_) < a 2 . 
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V. FINAL COMMENTS AND CONCLUSION 



In summary, we have obtained the approximate analytic relations for the relativistic 
energy spectra and the corresponding upper and lower spinor wave functions in the pres- 
ence of spherical scalar and vector reflectionless-type and RM potential models under the 
conditions of the spin and pspin symmetries. The resulting solutions of the wave func- 
tions are being expressed in terms of the generalized Jacobi polynomials or hypergeometric 
functions. Parametric generalization of the NU method is used. We have further used 
the recently introduced exponential approximation to deal with the spin-orbit centrifugal 
(pseudo-centrifugal) potential term. The most stringent interesting result is that the present 
spin (pspin) symmetric energy spectrum of the Dirac equation is noticed to be the same as 
the energy spectrum of the KG solution if V" = ±S (i.e., E = A = 0, C± = 0). We 
point out a possible remark of this result. The conditions that originate the spin and pspin 
symmetries in the Dirac equation are the same that produce equivalent energy spectra of 
relativistic spin-1/2 and spin-0 particles in the presence of spherical vector and scalar po- 
tentials. Obviously, the relativistic solution can be reduced to its non-relativistic limit by 
the choice of appropriate mapping transformations. Also, in case when spin-orbit quantum 
number k = ±1 (1 = 1 = 0), the problem can be easily reduced to the s-wave solution. We 
also find that when we let x — > —x, a — > ia, — > 1V2, the RM potential (1) turns into tRM 
potential (71) with real energy solution. We hope that, as in the nonrelativistic case (see, for 
example [58]), the relativistic model under consideration can be applied in molecular physics 
as well as nuclear physics. We stress that the present results should be useful in studying 
the rotation-vibration energy spectrum of low vibrational molecules of small amplitude and 
spin-orbit quantum nnmber k that are not large [51,58]. 
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FIG. 1: A plot of the Rosen- Morse-type potential (1) for three different cases V2 = 2V\, V2 = V\ 
and V 2 = Fi/3. 
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r (in units of 1/oc) 



FIG. 3: A plot of the Rosen-Morse potential (54) for three different sets of parameters (i) a = 5, 
a = 1, b = 1, (ii) a = 1, a = 5, b = 2, and (hi) a = 2, a = 10, b = 1. 
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FIG. 4: Plot of the trigonometric Rosen-Morse-type potential [see Eq. (71)] for various sets of 
parameters. 
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